Mult i walled carbon nanotube: Luttinger liquid or not? 
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We have measured IV curves of multiwalled carbon nanotubes using end contacts. At low voltages, 
the tunneling conductance obeys non-Ohmic power law, which is predicted both by the Luttinger 
liquid and the environment-quantum-fluctuation theories. However, at higher voltages we observe 
a crossover to Ohm's law with a Coulomb-blockade offset, which agrees with the environment- 
quantum-fluctuation theory, but cannot be explained by the Luttinger-liquid theory. From the 
high-voltage tunneling conductance we determine the transmission line parameters of the nanotubes. 
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Metallic carbon nanotubes are considered as outstand- 
ing realizations of strongly interacting, one-dimensional 
electron systems, i.e. Luttinger liquids (LL) |Q-[|. A 
Luttinger liquid is a paramagnetic metal without Fermi- 
liquid quasiparticles. Its basic charged excitations are 
plasmons which can be viewed as propagating electro- 
dynamic modes in a similar fashion as in any regular 
transmission line. Experimental evidence for LL behav- 
ior has recently been observed in single walled carbon 
nanotubes 0] as well as in multiwalled tubes (MWNT) 
H . The transmission line analogy, in turn, facilitates the 
connection of LL theory to the environment-quantum- 
fluctuation (EQF) theory fflM. This theory has been 
successful in explaining Coulomb blockade in normal tun- 
nel junctions Unlike the LL model, the EQF the- 
ory incorporates various factors, which makes it much 
more amenable to detailed experimental comparison, es- 
pecially in the case of resistive transmission lines. 

In this paper we present experimental results on 
the TV-curves for four metallic, arc-discharge-grown 
MWNTs. We analyze our results using both the EQF 
analysis HQ as well as the standard LL formulas [^|j9| . At 
small voltages both of these approaches predict for IV- 
curves a power law I cx V a+1 , which is also supported by 
the experiments yielding a + l = 1.23 ±0.1. At large volt- 
ages we find that only the EQF theory is applicable and 
obtain for the high frequency impedance Z = 1.3 — 7.7 
kf2. From the values of a and Z we can independently 
determine the kinetic inductance of the nanotubes and 
obtain consistent values of l^m = 0.1 — 4.2 nH/pm. 

A multiwalled nanotube consists of several concentric 
nanotubes. About one third of the tubes are expected to 
be metallic with quite large inter-layer capacitance. Ac- 
cording to our analysis, the large capacitance connects 
the inductive as well the resistive components of separate 
tubes in parallel. All metallic tubes take part in the con- 
duction at high frequencies, in contrast to the Aharonov- 
Bohm experiments of Ref. B where only the outermost 



layer contributed to the dc-resistance. The total electron 
density n of a MWNT is proportional to M, the number 
of metallic layers, and in each channel the Fermi veloc- 
ity of ID electron gas is vf = Trhn/AMm* = 8 • 10 5 m/s. 
Here m* is the effective mass of an electron, and we have 
taken into account the fact that each metallic layer has 
four independent ID conduction channels. 

A metallic nanotube, placed on a silicon substrate be- 
tween metallic contact lines, can be viewed as an inner 
conductor of a transmission line whose outer conductor 
is formed by nearby metallic bodies. The capacitance 
per unit length of the line is c = 27reeo/ln(r g /ro). Here 
e is the dielectric constant of the medium between the 
conductors, r is the outer radius of the nanotube, and 
r g is a distance from a metallic ground. The current 
carriers of the nanotube occupy ID conduction bands 
and, in contrast to the carriers in metallic wires, they 
have a low total density n, resulting in a large kinetic 
energy stored in the current flow. Therefore, the mag- 
netic inductance l m = iiq\d.(t g / r^) / 2~k , which is usually 
relevant for transmission lines, has to be replaced by the 
kinetic inductance l^m = m*/ne , since lki n 3> lm- In 
addition, for a ID plasmon in a nanotube the inverse 
compressibility dfi/dn = m*Vp jn of the neutral Fermi 
gas becomes comparable to the electrostatic inverse com- 
pressibility e 2 /c of the transmission line geometry. This 
can be taken into account by renormalization of the nan- 
otube capacitance into c: 
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and in the expression for the line impedance Z = V/I : 
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V is the electrochemical (not only electric) potential dif- 
ference, and Rk = h/e 2 is the quantum resistance. 

In the above classical electrodynamic analysis the ID 
plasmon modes are the only excitations of the nanotube 
transmission line. The LL model for an infinitely long 
MWNT [||,^| recovers the electrodynamic plasmon mode 
with v p i given by Eq. (Q) (cf. the expression after Eq. 
(4) in Ref. ||). But in addition to the plasmon mode, 
the Luttinger liquid has charge-neutral modes (a spin 
wave among them), which propagate with the velocity 
different from v p i and keep the total charge density con- 
stant. The Coulomb interaction, measured by the differ- 
ence of Vpi/vF from unity, suppresses the single electron 
density of state (DOS) p(E) — dn/dE near the Fermi 
level. The DOS is given by the Fourier component of 
the electron Green's function (ip(x,t)ip(x,0)^} (ip(x,t) is 
the electron operator) and is probed by the IV curve: 
dl/dV oc p{eV). At low energies p(E) cx E aL , where for 
an end-contacted infinitely long MWNT ||@ 

a L = {v p i/v F - 1)/4M . (4) 

In the limit of large M or no- interaction v p i/vf = 1, 
the TV-characteristics of a MWNT approach Ohm's law 
{oll = 0), i.e. the Luttinger liquid turns into a Fermi- 
liquid. 

Another approach, the EQF theory, considers the 
effects of environment quantum fluctuations on TV- 
characteristics under the conditions of Coulomb block- 
ade. A non-resistive, infinitely long nanotube acts as 
a dissipative environment, i.e., as a heat bath with 
which the tunneling electron can exchange energy Hf7j. 
The energy exchange is characterized by the function 
P[E), which is a Fourier component of the correlator 
(e'^^e"'^' '}, where <p(t) is the operator of the phase. 
At T = Q, P(E) is proportional to the second deriva- 
tive of the current: d 2 I/dV 2 oc P(eV). In the Coulomb 
blockade regime, i.e. when the voltage bias is less than 
e/Cx, where Ct is the capacitance of the tunnel con- 
tact, the EQF theory predicts that P(E) oc B"^ 1 
and I oc V aE+1 w ith a E = 2Re{Z}/R K . Using the 
impedance Z = \Jlkin/c of the nanotube, this yields 
the same power law as the LL theory in the large M 
limit. This fact, pointed out in Ref. fllPf , is not acciden- 
tal. In the LL picture the current is suppressed because 
there are no single electron quasiparticles, and the charge 
is transported by bosonic modes (plasmons). Although 
in a junction between 3D wires there are single-electron 
states available (in contrast to ID), a tunneling electron 
at V << e/Cx has not enough energy to get into them. 
As a result, the charge is transported again with ID plas- 
mons, which have similar properties for ID and thin 3D 
wires. 

On the other hand, one should expect a similarity be- 
tween p(E) and P{E), since both the operators, if)' and 
e~ l<p , which define these two functions, are creation op- 
erators for the charge e. But if the exponents ai and 



aE for the conductance coincide, the exponents ql and 
as — 1 for p(E) and P(E), respectively, differ by one. 
One can show |0, however, that similar relations con- 
nect the p{E)- and P(T?)-exponents with the impedance: 
a L = 2Rc{Z L }/R K - 1 and a E - 1 = 2Rc{Z}/R K - 1. 
But due to charge-neutral modes, which were not consid- 
ered in the EQF theory, the nanotube impedance differs 
in the LL theory from the impedance Z given by Eq. 
(§): Z L = Z + {AM - 1)R K /8M. The difference in the 
impedance compensates the difference in the relations 
connecting p(E) and P(E) with the conductance, and 
eventually in the large M limit both the theories predict 
the same exponent. 

But the two approaches differ in their predictions 
for high voltages. According to the EQF theory the 
power law is only valid in the Coulomb blockade regime 
ujCt « 1/Z. The relevant frequency u) ~ eV/fi in 
this inequality corresponds to the environment mode 
excited by a tunneling event; in our experiments this 
means frequencies up to about 20 THz. At high fre- 
quencies and voltages the environmental impedance Z is 
shunted by tunnel junction capacitance Ct and becomes 
(iuiCt + 1 Then the EQF theory gives the formula 
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This high-voltage asymptotics, characterized by the 
Coulomb offset e/2Cr and the "tail" voltage oc 1/V was 
experimentally studied and discussed by Wahlgren et al., 
and Penttila et al. within the horizon picture fll2] , |l3| |. 
In contrast to the EQF, in the LL approach the capaci- 
tance Ct of the tunneling contact is absent, and therefore 
this approach does not predict a crossover to the "tail" 
asymptotics given by Eq. (||). 

A summary of our four nanotube samples, each with a 
diameter of about 15 nm, is presented in Table I. For con- 
tact, we employed gold electrodes which were evaporated 
either prior to or after the deposition of nanotubes. De- 
position of nanotubes was done as described in Ref. . 
Mapping of nanotubes with respect to alignment marks 
as well as AFM micromanipulation was performed us- 
ing Park Scientific Instruments Autoprobe CP. Chrome 
or titanium (2-3 nm layer) was employed as an attach- 
ment layer before evaporating gold. Vacuum brazing at 
700 C for 30 sec was employed to lower the contact re- 
sistance in samples T1-T3. On the dilution refrigerator, 
the samples were mounted inside a tight copper enclosure 
and the measurement leads were filtered using 0.5 m of 
Thermocoax cable. 

Tunnel junction capacitances Ct — 31 — 111 aF and 
resistances Rt = 20 — 68 kfl (neglecting the tube re- 
sistance) were determined from asymptotic behavior by 
fitting Eq. (4) to the measured TV-curves. Owing to 
their relatively large size, the contacts to the nanotube 
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are not ideal and may cause a small uncertainty in the 
interpretation of the a values. Namely the EQF theory 
and the LL model in the strong interaction and M » 1 
limit predict two times smaller a values for the bulk than 
for the end contact 

Fig. 1 illustrates the low-voltage TV-curve of all four 
samples T1-T4. We are plotting the quantity V ff set — 
V — 1/-%* vs. V, which in the case of a power law yields 
a straight line with the slope a/(l + a). Only slight de- 
viation of linear behavior is seen at low voltages in Fig. 
1. This indicates that the Coulomb blockade of the is- 
land is rather weakly seen (except in T4). The linear 
behavior also implies that the two tunnel junctions be- 
come independent. At voltages \V\ > 5 mV, in spite of 
the additional wiggles, slight tendency toward saturation 
is observed in the data. This is consistent with the EQF 
picture, which predicts that at high voltages V ff se t must 
gradually approach to a Cr-dependent constant. By fit- 
ting a straight line through each data set at 4 < |V| < 7 
mV, we obtain a = 0.23 ± 0.1 (Table 1). Using Eqs. (2) 
and (3) we obtain lkin a = 0.2 — 1.7 nH//im for the kinetic 
inductance. Capacitance for the nanotube ~ 70 aF//im 
is estimated using the average of c = 2Ct/ '£ where we 
employ the total tube length C for scaling. 

Fig. 2 displays TV-curves measured at large voltages. 
In order to facilitate a direct comparison with the power- 
law dependence, we have plotted our results on a log-log 
scale. Our data are rather close to a single power law 
with small a but, at larger values of a (samples Tl and 
T3), there is a gradual approach toward a linear law as 
expected for a single junction in a resistive environment. 
Thus both figures give evidence that the environmental 
(EQF) theory is better suited for the analysis at high 
voltages. In fact, also the saturation observed by Bock- 
rath et. al. [|J can be explained by asymptotic approach 
toward Ohm's law. 

The plasma resonances, which one expects in finite 
nanotubes, are washed away in our samples. This gives 
a lower limit for the resistivity of the line, rC > Z/A. On 
the other hand, the LC-line model works over our volt- 
age range, i.e. r < Lolkim which results in an upper limit 
for r of the order of 1 k£!//im (at 1 mV). For compari- 
son, from the two-terminal resistance measurements we 
estimate that r < 20 kil/pm for our tubes at DC |l5| . 

One may argue that the poor agreement of the LL pic- 
ture with the high-voltage part of the IV curves could 
be reconciliated by including the junction capacitance 
Ct into the impedance, as is done in the EQF theory. 
There is, however, a conceptual problem to do it. The 
density of state p is expected to be a bulk property and, 
therefore, independent of Ct- Moreover, inclusion of Ct 
into the impedance, which determines p, does not help to 
match the LL picture with experimental results. The ca- 
pacitance Ct short circuits the environment impedance, 
and p(E) should decrease with E, like P{E) in the EQF 
theory. But since dl/dV oc p(E), in contrast the EQF 



theory where dI 2 /dV 2 oc P(E), this yields a high- voltage 
plateau (voltage-independent current), but not Ohm's 
law with Coulomb offset. Introduction of a proper high- 
energy cut-off in the LL model could explain a cross-over 
to Ohm's law, but not a Coulomb offset. We expect this 
cut-off to be larger than the region of our analysis which 
is bounded by the presence of higher transverse modes 
above 50 mV. 

Fits, based on Eq. (JsJ) , fall on top of the experimental 
data in Fig. 2. In the fitting, we assume that the junc- 
tions at the ends of the tube are symmetric and, in fact, 
T vs. V/2 is fitted to the single junction formula. We 
also tried to incorporate a cubic background, r)V 3 in the 
fitting, which was found essential in Al-samples because 
of the deformation of the tunnel barrier at high voltages 
|l3f . Surprisingly, the cubic term was found negligible in 
all our nanotube samples. Our fits yield a characteristic 
impedance of Z = 1.3 — 7.7 kSl for the resistive environ- 
ment. These results depend slightly on the measurement 
polarity (see Table I). Finally, using Eqs. (2) and (3), we 
obtain for the kinetic inductance l^m = 0.1 — 4.2 nH//im. 

Table I contains parameters obtained both from the 
EQF analysis for a LC-transmission line as well as from 
the power law exponents according to the LL model. The 
results of the two methods overlap each other; the scat- 
ter of the power law analysis is slightly smaller than 
that of the environmental analysis. In addition, we 
checked that the temperature dependence of the mea- 
sured conductance, dl/dV oc T a , yielded consistent val- 
ues of a — 0.25 ±0.1. As a final result of all our 
determinations we quote the median value Ikin = 0.5 
nH//im. If we compare this with the theoretical pre- 
diction Ikin = Rk /8Mvp, we conclude that the average 
number of conducting layers in our nanotubes is 8 and 
the large variation of the inductance may come from the 
variation in M. The average value of 8 indicates that 
about every 3rd layer in our nanotubes is metallic. 

To conclude, on the basis of our experimental results 
we argue that, at high voltages, the environmental the- 
ory gives a better account of transport measurements of 
multiwalled nanotubes than the Luttinger liquid picture, 
because the tunnel junction capacitance is neglected in 
the Luttinger liquid theory. At lower voltages, no dis- 
tinction between these two theories can be made. Due to 
their large kinetic inductance, nanotubes provide an ex- 
cellent high-impedance environment for normal junctions 
at high frequencies, which is crucial for single-electronics 
phenomena. As the kinetic inductances of different layers 
are in parallel in MWNT, these phenomena will be more 
pronounced in single walled carbon nanotubes. 
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TABLE I. Summary of our samples of single, metallic 
tubes T1-T4. C1/C2 denotes the length of the tube over the 
metallic leads/the length of the free standing section that is 
hanging ~ 20 nm above the substrate. For sample T4, £2 
specifies the length in contact with Si02- In the resistiv- 
ity ratio rr = i? (4.2K)/i? (290K), the values of R have 
been obtained from slopes of TV-curves at I = 0. T denotes 
the measurement temperature for the data in Figs. 1 and 2. 
The junction capacitance Ct and tunneling resistance Rt are 
taken from the IV-cwcve fits using Eq. (5). These fits also 
yield the lumped-element environmental impedance Z which 
is slightly different for positive and negative voltages. Kinetic 
inductance obtained from Z using c = 70 aF/fim and Eq. (3) 
is given by lkin z ■ The Luttinger liquid power a is determined 
at 4 < [V| < 7 mV. The kinetic inductance lkin a is an esti- 
mate obtained using Eq. (3). Voltage range of the tail fits is 
given in the last column. 



C1/C2 rr T R T Ct 
fim K kO aF 

Tl 0.7/0.5 1.7 4.2 25 31 
T2 0.3/0.8 3.0 0.1 20 33 
T3 0.8/0.6 4.6 4.2 46 37 



Z lkin z a 

kfl mH/m 

3.5/7.7 0.9/4.2 0.30 

1.3/2.3 0.1/0.4 0.12 

2.0/4.8 0.3/1.6 0.32 



lkin a Range 
mH/m mV 

1.1 10-50 
0.2 15-50 

1.2 10-50 



T4 1.4/2.3 3.0 0.1 68 111 1.8/2.7 0.2/0.5 0.17 0.3 7-20 
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FIG. 1. Offset voltage V a f fset = V - I / vs. V for all 
our samples T1-T4; the power law behavior / oc V a+1 yields a 
straight line in this kind of a plot. The dashed lines illustrate 
linear fits made in the range 4 < \V\ < 7 mV. The effect of 
Coulomb blockade near zero is seen to be small except for the 
sample T4. 
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FIG. 2. High voltage IV-curves (both positive and negative 
polarities) on a log-log plot. The dashed line illustrates linear 
behavior (a = 0). For details, see text. 
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